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X-RAY TIMING OF THE ENIGMATIC NEUTRON STAR RX J0720.4-3125 



(N 
O 
O 
(N 

Or 

00 

(N 
> 

(N 
O 

in 
o 

(N 
O 



D. L. Kaplan, S. R. Kulkarni 

Department of Astronomy, 105-24 California Institute of Technology, Pasadena, CA 91125, USA 
dlk@astro.caltech.edu, srk@astro.caltech.edu 

M. H. VAN Kerkwijk 
Sterrenkundig Instituut, Universiteit Utrecht, Postbus 80000, 3508 TA Utrecht, The Netherlands 

M.H.vanKerkwijk@phys.uu.nl 
AND 

H. L. Marshall 

Center for Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 

hermanm@space.mit.edu 
Accepted for publication in ApJL 

ABSTRACT 

RX J0720.4— 3125 is the third brightest neutron star in the soft X-ray sky and has been a source of 
mystery since its discovery, as its long 8-s period separates it from the population of typical radio pulsars. 
Three models were proposed for this source: a neutron star accreting from the interstellar medium, an 
off-beam radio pulsar, or an old, cooling magnetar. Using data from Chandra, ROSAT, and BeppoSAX 
we are able to place an upper limit to the period derivative, \P\ < 3.6 x 10~ 13 s s" 1 (3-cr). While 
our upper limit on P allows for the accretion model, this model is increasingly untenable for another 
similar but better studied neutron star, RX J1856.5— 3754, and we therefore consider the accretion 
model unlikely for RX J0720.4-3125. We constrain the initial magnetic field of RX J0720.4-3125 to be 
< 10 14 G based on cooling models, suggesting that it is not and never was a magnetar, but is instead 
middle-aged neutron star. We propose that it is either a long-period high-magnetic field pulsar with 
P ~ 1CP 13 s s" 1 similar to PSR J1814— 1744, or a neutron star born with an initial period of w 8.3 s 
and P ~ 10 -15 s s _1 . The proximity of RX J0720.4— 3125 is strongly suggestive of a large population 
of such objects; if so, radio pulsar surveys must have missed many of these sources. 

Subject headings: pulsars: individual (RX J0720.4— 3125) — stars: neutron — X-rays: stars 
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1. INTRODUCTION 

RX J0720.4-3125 was discovered by Haberl et al. (1997) 
as a soft (kT ~ 80 eV), bright X-ray source in the ROSAT 
All-Sky Survey. Given its very low hydrogen column 
density (Njj ~ 1 x 10 20 cm -2 ), nearly sinusoidal 8.39-s 
pulsations, relatively constant X-ray flux, and very faint 
(B = 26.6 mag), blue optical counterpart (Kulkarni & 
van Kerkwijk 1998; Motch & Haberl 1998), it was classi- 
fied as a nearby, isolated neutron star. 

As one of the closest (d ~ 300 pc; Kaplan, van Kerkwijk, 
& Anderson 2002) neutron stars, RX J0720.4-3125 occu- 
pies a central position in our study of these objects. How- 
ever, the long period is puzzling, and has led to three mod- 
els: an old, weakly magnetized neutron star accreting mat- 
ter from the interstellar medium (Wang 1997; Konenkov 
& Popov 1997); a middle-aged pulsar with ~ 10 12 G mag- 
netic field whose radio beams are directed away from the 
Earth (Kulkarni & van Kerkwijk 1998); or an old mag- 
netar (neutron star with magnetic field > 10 14 G; Dun- 
can & Thompson 1992) that is kept warm by the decay 
of its strong magnetic field (Heyl & Hernquist 1998; Heyl 
& Kulkarni 1998). These models predict different period 
derivatives: P < 5 x 10~ 15 s s" 1 , P ~ 10 -15 -10" 13 s s" 1 , 
and P > few x 10" 13 s s _1 , respectively 

Motivated thus, we undertook timing observations of 
RX J0720.4— 3125 using the Chandra X-ray Observatory 
(CXO), supplemented with analysis of archival data from 
ROSAT and BeppoSAX. 

After submission of this paper, we became aware of the 
work of Zane et al. (2002) reporting a timing analysis of 



RX J0720. 4-3125. A "Notes added in manuscript" section 
regarding that analysis can be found at the end of the 
manuscript. 

2. OBSERVATIONS 

The primary data consist of two sets of observations ob- 
tained from Chandra: one using the HRC in the spectro- 
scopic mode (HRC-S) with the Low Energy Transmission 
Grating (LETG), and one using ACIS in the continuous 
clocking (CC) mode. The primary and archival datasets 
are summarized in Table 1. 

We processed the HRC-S data using the standard 
pipeline 1 and extracted th order events from a circle with 
radius 10 pixels (l'/3). For the ±l st orders, we extracted 
events from a region 0.0006° wide in the cross-dispersion 
direction (the tg_d coordinate) and from 0.08° to 0.35° 
along the dispersion direction (the tg_r coordinate). We 
extracted events from the ACIS data within ±1" of the 
source. We then used the axBary program to barycenter 
the events in both these datasets. 

The best fit position for RX J0720.4-3125, found by av- 
eraging the th order data from the three HRC-S datasets, 
is (J2000) a = 07 h 20 m 24?96, S = -31°25'49'/6, with rms 
uncertainty of « 0'.'6 in each coordinate due to CXO as- 
pect uncertainties. This is consistent (l'/4 away) with the 
optical position (Kulkarni & van Kerkwijk 1998). The X- 
ray source appears unresolved and its profile is consistent 
with that of a point source (half-power radius of « 0'.'5). 

For the ROSAT HRI data, we extracted the events 

http:/ /asc. harvard.edu/ciao/threads/spectra_letghrcs/ 
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Table 1 
Summary of Observations 



Date 


MJD 


Exp. 


Span 


Counts 


BG 


Facility 


Instrument/ 


Period 


TOA - MJD 50000 d 




(day) 


(kscc) 


(kscc) 




Counts' 1 




Mode b 


(s) 


(TDB day) 


1993-Sop-27 


49257.2 


3.2 


12.0 


5800 


22.8 


ROSAT 


PSPC 


8.3914(4) 


-742.745297(3) 


1996-Nov-03 


50390.9 


33.7 


65.7 


12662 


79.0 


ROSAT 


HRI 


8.39113(6) 


391.300750(2) 


1997-Mar-16 


50523.1 


18.1 


99.4 


407 


15.4 


BeppoSAX 


LECS 


8.39103(9) 


523.705635(4) 


1998-Apr-20 


50923.2 


8.1 


460.3 


3074 


17.1 


ROSAT 


HRI 


8.391114(14) 
8.39111(2)° 


925.688213(5) 


2000-Feb-01 


51575.3 


5.4 


305.5 


929 


1.3 


CXO 


HRC-S+LETG 


1577.039569(2) 










671 


127.0 


CXO 


HRC-S+LETG ±1 


2000-Fob-02 


51576.1 


26.3 




4584 
3027 


5.2 
454.0 


CXO 
CXO 


HRC-S+LETG 
HRC-S+LETG ±1 






2000-Feb-04 


51578.7 


6.1 




1119 
687 


1.2 
119.5 


CXO 
CXO 


HRC-S+LETG 
HRC-S+LETG ±1 






2001-Doc-04 


52247.7 


15.0 


168.6 


31746 


229.8 


CXO 


ACIS-S3/CC 


8.391119(12)° 


2248.6768200(8) 


2001-Dec-05 


52248.2 


10.6 




22825 


155.8 


CXO 


ACIS-S3/CC 


2001-Dec-06 


52249.6 


4.1 




8786 


61.4 


CXO 


ACIS-S3/CC 







^Background counts scaled to the source extraction area. b HRC-S+LETG indicates order 0; HRC-S+LETG ±1 indicates orders ±1. 
c Valucs in parentheses arc 1-cr errors in the last decimal digit. d TOA is defined as the maximum of the folded lightcurvc nearest the middle 
of the observation, as determined from the best-fit sine wave. The ACIS/CC times were corrected for spacecraft motion following Zavlin et al. 
(2000). °A11 pointings for each of the Chandra HRC-S and Chandra ACIS datasets were processed together. 



within a circle of radius 45 pixels (22'.'5) centered on the 
source. We used a circle of radius 200 pixels (100") for 
the PSPC data. These events were barycentered using the 
f tools programs abc and bet and corrected to Barycen- 
tric Dynamical Time (TDB) according to Cox (2000, p. 
14). 

We extracted the BeppoSAX LECS events within a cir- 
cle with radius of 25 pixels (200") centered on the source 
and restricted to those with pulse-invariant (PI) ampli- 
tudes that were less than 90 (energies < 0.95 kcV), in order 
to maximize the signal-to-noise. Finally we barycentered 
the events with the SAXDAS tool baryconv. 

3. TIMING ANALYSIS 

For each dataset, we computed Z\ power spectra around 
the known 8.39-s period. Specifically, we explored the pe- 
riod range from 8.376 s to 8.405 s in steps of 7 /xs (oversam- 
pling by factors of 20-800 relative the nominal step-size 
of Pg/AT, where Po = 8.39 s is the approximate period 
and AT is the span of the dataset from Table 1). As can 
be seen from Figure 1 all but the Chandra HRC-S and 
ROSAT HRI-2 datasets yielded unambiguous period esti- 
mates. For the HRC-S and HRI-2 sets the period estimates 
are ambiguous because the large gaps in the observations 
result in strong side-lobes. In Figure 2, we display the best- 
fit periods for the unambiguous determinations as well as 
possible periods for the HRC-S and HRI-2 datasets. 

As can be seen from Figure 2, the ambiguity of the HRC- 
S and HRI-2 datasets can be resolved provided we assume 
(reasonably) that the period evolves smoothly with time. 
Our choice of period (for HRC-S and HRI-2) and the best 
fit periods (for the other datasets) are shown in Table 1. 
The errors on the periods were determined using the an- 
alytical expression from Ransom (2001). While that ex- 
pression was derived for FFT power spectra, Z\ power 
spectra have the same statistics (both are exponentially 
distributed) so the same relations should apply (we have 
verified this with numerical simulations). We also show in 
Table 1 times-of-arrival (TO As) for each of the datasets. 

The data in Table 1 are consistent with there being 
no measurable P: fitting for a linear spin-down gives 
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Fig. 1. — Z'l pcriodograms for the datasets listed in Table 1. 
For each dataset, the Z\ power is normalized so as to have unit 
mean (when no signal is present). Given that the statistics of Z\ 
are exponential it follows that the variance is also unity. Note the 
different vertical scales which reflect the differing significance levels 
of the detections. 

P = 8.391115(8) s at MJD 51633 and P = (1 ± 12) x 
10~ 14 s s _1 , with x 2 — 1-54 for 4 degrees-of-freedom 
(DOF). If, instead, we fit only for a constant period, we 
find P = 8.391115(8) s, with X 2 = 1-64 for 5 DOF. There- 
fore we can constrain the secular period derivative to be 
\P\ < 3.6 x 10~ 13 s s- 1 (3-cr). 

The folded lightcurve is largely sinusoidal, with an rms 
pulsed-fraction (the rms of the lightcurve divided by the 
mean) of 8% for both Chandra datasets. However, this 
pulsed-fraction is energy-dependent: the fraction rises 
with decreasing energy (see Figure 3), in agreement with 
the XMM analysis (Paerels et al. 2001). 

4. DISCUSSION 

Our upper limit of \P\ < 3.6 x 10~ 13 s s _1 is sufficiently 
high that we cannot meaningfully constrain the accretion 
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model, for which we expect P < 5 x 10 15 s s 1 (the limit 
is the case where all of the required M of ~ 10 12 g s _1 
couples to the neutron star at the corotation radius, giv- 
ing maximum torque per unit mass). We note that the 
accretion model is no longer viable for another similar 
but better studied isolated neutron star, RX J1856.5— 3754 
(van Kerkwijk & Kulkarni 2001). Regardless, the accretion 
model is best confronted by measuring the proper motion 
and distance, and looking for evidence of sufficiently dense 
ambient gas (deep Keck Ha imaging and HST astrometric 
observations are in progress). 
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Fig. 2. — Period measurements for RX J0720. 4-3125, using 
the data from Figure 1. As explained in the main text the Chan- 
dra HRC-S and ROSAT HRI-2 measurements are ambiguous ow- 
ing to large gaps in the data. Probable periods arc displayed. 
The best-fit constant period model is shown by the dashed line: 
P = 8.391115(8) s. 

However, we can constrain the pulsar and magnetar 
models. We can draw four inferences common to both 
models. First, the spin-down luminosity, E = Ilulu < 
2.4 x 10 31 erg s _1 ; here, u = 2-k/P. Second, in the 
framework of a simple (vacuum magnetic dipole radia- 
tion) pulsar model, the physical age is roughly approxi- 
mated (provided the current spin period is much larger 
than that at birth and that the magnetic field does not 
decay significantly) by the so-called characteristic age: 
r c = P/(2P) > 4xl0 5 yr. Third, the strength of the dipole 
field is B = 3.2 x 10 19 (PP) 1 / 2 < 6 x 10 13 G. Fourth, we 
assume that the X-ray emission (well described by a black- 
body; Haberl et al. 1997; Paerels et al. 2001) is cooling 
flux from the surface. The bolometric cooling luminosity 
is L cool :=y 2 x 10 32 d| 00 ergs s _1 (Haberl et al. 1997), using 
the distance estimate of 300d 30 o pc derived by scaling from 
RX J1856.5-3754 (Kaplan et al. 2002). 

Knowledge of i coo i enables us to estimate the cooling 
age, icooi, of RX J0720.4— 3125. Magnetic fields, especially 
intense B fields such as those proposed for magnetars, can 
profoundly influence the cooling of neutron stars. To this 
end, we use the curves of L vs. t from Heyl & Kulkarni 
(1998) and find t cool w (5 - 10) x 10 5 yr, assuming a 50% 
uncertainty in the distance and with only a slight depen- 
dence on B. This age is consistent with the characteristic 
age derived above. 
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Fig. 3. — RMS pulsed-fraction (see § 3 for discussion) for different 
energy bins, from the Chandra ACIS data. The overall pulsed frac- 
tion is 8.1 ± 0.4%. Each bin was chosen to have the same number 
of total counts. 



In the magnetar model (we assume that the £?-ficld de- 
cay is dominated by the slower irrotational mode; see Heyl 
& Kulkarni 1998) the expected B field at about 10 6 yr is, 
w 2 x 10 14 G, well above the upper limit obtained from 
our P limit. Models that are consistent with both L coo i 
and our limit on P (and thus an upper limit on the current 
value of B) are those with initial B < 10 14 G. Based on 
this, we conclude that RX J0720.4— 3125 is not a magne- 
tar, motivating us to consider the pulsar model. 

Earlier, Kulkarni & van Kerkwijk (1998) did not ac- 
cept the radio pulsar model because in 1998 there were no 
radio pulsars with such long periods. However, over the 
past four years we have come to appreciate the existence of 
pulsars with B > 10 13 G (Camilo et al. 2000; see also Fig- 
ure 4). In particular, the parameters of RX J0720.4— 3125 
are not too dissimilar to those of PSR J1814— 1744, which 
has P w 4 s and P w 7.4 x 10~ 13 s s _1 . Thus, the past ob- 
jections against the radio pulsar model are no longer ten- 
able, and RX J0720.4— 3125 seems fully compatible with 
being an off-beam high-P pulsar. If that is the case, then 
we expect P ~ 10~ 13 s s _1 , a value that we should be able 
to measure in the near future. 

A separate possibility is that RX J0720. 4-3125 is 
an off-beam pulsar with age compatible with t coo \, but 
with a conventional (~ 10 12 G) magnetic field and P ~ 
10~ 15 s s- 1 (and therefore r c - 10 8 yr). RX J0720. 4-3125 
could then be similar to the 8.5-s, 2 x 10 12 -G pulsar 
PSR J2144-3944 (Young, Manchester, & Johnston 1999; 
see Figure 4). With a braking index of 3, the age of a 
pulsar is t — t c (l — (Pq/P) 2 ^, where Po is the initial 
spin-period. If RX J0720.4-3125 does have T ~ icooi and 
B ~ 10 12 G, we find and P « 8.3 s, very close to P. 
Such as a pulsar would be an example of the "injection" 
hypothesis (Vivekanand & Narayan 1981), where pulsars 
are born with initial spin periods Po 3> 10 ms (as for the 
Crab). Such long initial periods are allowed and perhaps 
expected in some models of neutron-star formation (e.g., 
Spruit & Phinney 1998), where the precise initial period 
depends very sensitively on the details of the formation 
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mechanism and may range over four orders of magnitude. 
While there are a few pulsars whose characteristic ages are 
factors of 10-100 times the ages derived from supernova 
remnant associations (Pavlov et al. 2002), this would be 
the first case for a source with P > 1 s. 
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Fig. 4. — P-P diagram, showing only P > 1 s and P > 
10 — 16 s s — 1 . Radio pulsars are plotted as points, magnetars as 
diamonds. RX J0720.4— 3125 is an upper limit. The magnetar 
IE 2259+58.6 is circled, as are the high-B pulsars PSR J1830-1135 
and PSR J1814-1744, and the long-period pulsar PSR J2144-3944 
(Young et al. 1999). A version of the so-called "death line" is marked 
by the dotted line. The sloping solid lines are lines of constant dipolc 
magnetic field B dipo i = 3.2 X 10 19 (PP) 1 / 2 G, while the dashed lines 
are those of constant characteristic age t c = P/(2P). 

We make the following parenthetical observation: for 
most known pulsars the X-ray pulsed-fraction (largely) 
increases with photon energy (Perna, Heyl, & Hcrnquist 
2001), whereas for RX J0720.4-3125 we sec the opposite 
effect. However, for the pulsars (e.g., PSR B0656+14, 
PSR B1055— 52), the X-ray luminosity has a strong, 
highly-pulsed, non-thermal component with ^x,non-th ~ 
lO -3 ^ (Becker & Trimmer 1997). Furthermore, in such 
objects heating of the polar caps by pulsar activity (prob- 
ably dependent on E) is likely significant. The interplay 
of these components with the viewing geometry can re- 
sult in the large range of observed phenomena (e.g., Perna 
et al. 2001). For RX J0720.4-3125, though, with its small 
E there is little reason to expect a strong non-thermal 
contribution or a hot polar cap (although there must be 
some inhomogeneities to give the observed pulsations). 
We conjecture that the increase in the pulse fraction with 
decreasing photon energy is primarily due to the absence 
the additional components. 

5. CONCLUSIONS 

In this Letter, based on X-ray timing data and cool- 
ing models, we argue that the nearby soft X-ray source 
RX J0720.4— 3125 is not a middle-aged magnetar but is 
likely a 10 6 -yr off-beam pulsar. To accommodate its age 
and long period we speculate that it either has B > 10 13 G 
or was born with P ~ 8.3 s, a very surprising result as 
both source types are, at present, considered to be rare. 



We now consider the larger ramifications of our conclu- 
sions. 

A volume-limited sample of neutron stars offers us an 
opportunity to sample the diversity of such sources. In 
this respect, soft X-ray surveys provide the best such sam- 
ples since all neutron stars — normal radio pulsars, high-B 
pulsars, magnetars, and the mysterious Cas-A-like neu- 
tron stars — will cool through soft X-ray emission well 
into their middle ages. Indeed, this expectation is borne 
out by the local sample: pulsars such as PSR B0656+14, 
Geminga 2 , RX J0720.4— 3125, a youngish magnetar (see 
below), and finally the very mysterious RX J1856.5— 3754, 
of which we know nothing other than it is a cooling neu- 
tron star 3 . 

The number of neutron stars belonging to a given class 
depends not only on the sensitivity of the X-ray sur- 
vey but also on the product of the birth rate and the 
cooling age. Thus, for example, magnetars with their 
longer-lasting cooling radiation may dominate the local 
population despite a lower birthrate (Heyl & Kulkarni 
1998). This bias and the long period led us to specu- 
late that RX J0720.4— 3125 was an old magnetar, a con- 
clusion we have now refuted. In contrast, the soft ther- 
mal X-ray source RX J1308. 8+2127, with P = 5.2 s and 
P ~ 10~ n s s _1 (Hambaryan et al. 2002) appears to have 
a magnetar-strength field. 

The proximity of RX J0720.4— 3125 argues for a sub- 
stantial Galactic population of similar sources, but very 
few such radio pulsars are known. The cause of this 
paucity is that radio surveys select against long-period pul- 
sars, especially those with B > 10 13 G, in several ways. 
(1) The beaming factor is known to decrease with increas- 
ing period, reaching 3% at P ~ 10 s (Tauris & Manchester 

1998) . (2) As can be seen from Figure 4, the lifetime of 
a radio pulsar decreases with increasing B: a B ~ 10 12 G 
neutron star crosses the the radio death line at ~ 10 8 yr 
whereas a B ~ 10 13 G pulsar dies at ~ 2 x 10 7 yr. The loss 
of throughput of a pulsar survey for a 5-s pulsar relative to 
a 1-s pulsar from these two effects alone is nearly one order 
of magnitude. (3) The true loss is even greater since long- 
period signals are frequently classified as interference (we 
note that population models do not constrain the popula- 
tion of long-period pulsars [Hartman et al. 1997], mainly 
due to reasons 1 and 2). Young high-B pulsars in super- 
nova remnants would almost certainly create visible ple- 
rion nebulae due to their high E's, while long-period in- 
jected pulsars of similar ages would be invisible except for 
their cooling radiation (without assuming that the radio 
beams are directed toward the Earth). Thus, injected pul- 
sars detectable only via X-ray emission may be present in 
many "hollow" supernova remnants (i.e. those without vis- 
ible plerions). Radio pulsar searches better tuned to long 
periods and very deep radio and X-ray searches for young 
pulsars in supernova remnants may uncover the postulated 
class of long-period sources. 

We thank M. Cropper and S. Zane for discussions on 
the TOAs and for alerting us to the different time sys- 

2 Presumably a standard pulsar that is not beamed toward us. 

' 5 It is further worth noting that the sample of soft X-ray neutron 
stars has at least three long-period objects (Haberl, Pietsch, & Motch 

1999) . 
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terns used by different programs. We thank D. Frail and 
J. Heyl for useful discussions, and we thank an anony- 
mous referee for constructive comments. We have used 
the NASA-maintained HEASARC web site for archival 
data retrieval and subsequent analysis. D.L.K. holds a 
fellowship from the Fannie and John Hertz Foundation, 
and his and S.R.K.'s research are supported by NSF and 
NASA. M.H.v.K. is supported by a fellowship from the 
Royal Netherlands Academy of Arts and Sciences. 

Notes Added In Manuscript After we submitted our pa- 
per we became aware of a paper by Zane et al. (2002) 
reporting timing analysis of RX J0720.4— 3125. Our pe- 
riod determinations of the archival data (PSPC, HRF1, 
LECS, HRF2, HRC-S) are in excellent agreement with 
those of Zane et al. Both papers also report new de- 
terminations, which are: Chandra ACIS-S3 (our paper) 
and XMMa (2000 May 13) and XMMb (2000 November 
21; both from Zane et al.). We restricted our analysis to 
an incoherent combination of the various datasets, i.e. we 
looked for secular evolution of the period determined from 
each observation separately. We did not attempt to phase 
connect the datasets. Zane et al. do present a coherent 
analysis, using the archival and XMM data. However, in 
our opinion such an analysis is premature and not robust. 
First, it is premature, because the XMM derived periods of 
known pulsars have fractional errors \AP/P\ ranging from 
1.9 x 10~ 7 to 1.2 x 10~ 5 (as reported by the XMM calibra- 
tion team; Kuster et al. 2002). This error alone may result 
in systematic uncertainties as high as P of 6 x 10~ 12 over 
the 6-month duration of the XMM datasets and 5 x 10~ 13 s 
s^ 1 over the entire span of the observations. Second, it is 
not robust, as phase connection without any ambiguity 
requires that the datasets be separated by time intervals 
less than "coherence" timescale, ~ P 2 /(Jp (where op is 
the uncertainty in the measurement of P) and none of the 
datasets (including the XMM datasets) satisfy this condi- 
tion. We note that neither of the two primary solutions 
from Zane et al. fits the TO As in Table 1. 
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